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Abstract 

We propose a protocol of quantum energy teleportation that trans- 
ports energy in spin chains to distant sites by only local operations 
and classical communication. The protocol uses ground-state entan- 
glement and localized negative-energy excitation, and the energy is 
teleported without breaking any physical laws including causality and 
local energy conservation. 
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1 Introduction 



Quantum entanglement opens the door to several interesting quantum 
phenomena, including quantum teleportation (QT) pQ, by which an unknown 
quantum state can be teleported to a distant place by local operations and 
classical communication (LOCC). However, QT cannot transport the excita- 
tion energy of the transported states. For example, let us imagine that Alice 
sends to Bob a spin-up state of a spin in an external uniform magnetic field 
parallel to the z axis. The Hamiltonian is given by 

H s = go z 

with a positive constant g. The ground state of the spin is the spin-down 
state with energy —g and the energy of the excited spin-up state is +g. A 
contracted state at Alice or Bob of a Bell pair shared by Alice and Bob for QT 
is the maximal entropy state, which has zero energy on average. Therefore, 
both Alice and Bob must first supply energy +g to their spins on average to 
create the Bell pair from two spins in the ground state by a global operation. 
Bob has to supply an additional +g energy to his spin to receive the spin-up 
state from Alice's site by QT. Hence, the total energy input at Bob's site is 
given by +2g, which should be locally prepared by Bob. 

Does the above- described protocol of QT imply that energy cannot be 
transported by LOCC? Amazingly, the answer is "no" . Quantum mechanics 
allows energy transport by LOCC. In this paper, we propose a protocol of 
quantum energy teleportation (QET) in spin chains which transports energy 
only by LOCC, using local excitations with negative energy and ground- 
state entanglement, thereby respecting fundamental physical laws including 
causality and local energy conservation. 

Spin chains are composed of many spins arrayed in one dimension. Short- 
range interactions exist between the spins, and the Hamiltonian is given by 
a sum of these local interaction terms. Due to the interactions, complicated 
entanglement among the spins can occur even in the ground state. Spin 
chains have recently attracted much attention in the context of quantum 
information theory because they can be applied to the short transmission of 
quantum states [2J. Spin-chain entanglement is also helpful for investigating 
the complicated physical properties of the ground state of spin chains [3]. 
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Negative localized energy plays an essential role in the protocol. In this 
paper, we define the zero values of energy density by the expectation values 
of the ground state. We then consider why there are negative-energy- density 
regions in the spin chains. In quantum physics, there remain local quantum 
fluctuations even in the ground states of spin chains. By linearly superpos- 
ing the eigenstates of the total Hamiltonian, we can suppress the quantum 
fluctuations more strongly in a local region via quantum interference, as 
compared with the ground state. Taking a definition such that the energy 
density in the ground state is zero, the energy density becomes negative 
in a region where the quantum fluctuation is more greatly suppressed. We 
might be concerned that there are states with energy values lower than the 
ground state. However, it should be stressed that, even though some regions 
have negative energy density, the total energy of the system is always non- 
negative. Hence, there exist no states with energy lower than the ground 
states. Negative energy effects in relativistic field theory have long been 
investigated [4J. However, they have not often been applied to condensed 
matter physics, quantum optics or quantum communication. As an excep- 
tional example, a fundamental lower bound of actuating energy for photon 
switching has recently been derived for input waves with negative energy 
density by a gedanken experiment [5]. 

In this paper, we present a protocol for near-critical two-level spin chains 
with nondegenerate ground states and large correlation lengths. In the anal- 
ysis we assume that the number of spins is quite large, but finite. However, 
the results obtained in this paper can also be applied to the case of infinite 
degrees of freedom as long as an infinite limit of the number of spins is justi- 
fied. We also concentrate on short time scales in which dynamical evolution 
induced by the Hamiltonian H is negligible. We denote the difference be- 
tween the largest eigenvalue and the smallest eigenvalue of H by AE. The 
timescale t considered is assumed to satisfy 

(<< ^' « 

Assuming this condition, it is valid to treat the time evolution operator as 
exp [—itH] ~ /. It should also be noted that the condition in eq. (JTJ can be 
weakened if a finite amount of energy Ei n less than AE is input to the spin 
chain by the energy sender as follows: 

* « ^- (2) 



2 



On the other hand, we also assume that LOCC for the spins can be repeated 
many times even in a short time interval. Taking the site number difference 
between the two parties in the protocol as An and the lattice spacing between 
nearest-neighbor sites as a, the time scale condition for many-round LOCC 
is expressed as 



where c is the light velocity. By taking the nonrelativistic limit c — > oo, 
the relation in eq. ([3]) always holds. The proposed protocol is presented for 
one-dimensional qubit chain systems. Extending the protocol to spin chain 
systems with larger spins and larger dimensions is a straightforward task. 

The paper is organized as follows. In section 2, we show that negative 
localized energy density naturally appears in the systems under study. In 
section 3, local energy conservation of the spin chain systems is discussed. 
In section 4, we discuss the measurement of local observables for the ground 
state, which is an important consideration for the QET protocol. In section 
5, we propose the QET protocol. We conclude the paper in section 6. 



2 Ground- State Entanglement and Negative 
Energy Density 



Let us consider a spin chain with a nondegenerate ground state. The 
Hamiltonian is given by a sum of semi- local components T n : 



t ^> aAn/c, 



(3) 




(4) 



n 



Here T n are Hermitian operators given by 



n+L 




(5) 



7 m=n—L 
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where Om is a local Hermitian operator at site m and the integer L denotes 
the interaction range. If we take L — 1, the nearest neighbor interaction can 
be treated. For example, the interaction of the Ising model with a transverse 
magnetic field satisfies L = 1 and has T n such that 

T n = ~ba z n - \ol « +1 + <_ x ) - e, (6) 

where b, h and e are real constants. T n of the example has three terms in the 
right-hand side in eq. ([5]) with 7 = 1,2, 3. These operators are given by 



w n-l 

^ n 

n (n,l) 
w n+l 

w n-l 


= -K- 
= i, 
= /, 




h X 
= -Ifr, 


n M 

w n+l 
n (n,3) 


= a n+l) 
= °n-li 


r)(n,3) 


h x 

= ~^n, 


n (n,3) 


= I. 



T n describes the local energy density at site n. The ground state \g) is an 
eigenstate for the lowest eigenvalue Eq of H. However, it is not guaranteed 
in general that \g) is an eigenstate of T n . Without changing the dynamics at 
all, it is always possible to shift T n by adding constants. Hence, we are able 
to redefine T n to satisfy 

(g\T n \g) = o. (7) 

Then, the shifted value of Eq becomes zero as follows: 

E = (g\H\g) = J2(9\Tn\g)=0. 

n 

From the above, it is satisfied that 

H\g) = (8) 



4 



without any loss of generality. By this redefinition, H becomes a non-negative 
operator and satisfies 

Tr [pH] > (9) 

for an arbitrary quantum state p. This choice of the energy origin simplifies 
the spin chain analysis, because we will often consider the energy difference 
between an excited state and the ground state. 

If each T n is an exact local operator at site n, all T n can be simultaneously 
diagonalized and the ground state \g) is separable and an eigenstate for the 
lowest eigenvalue of each T n . In such a situation, T n becomes non-negative. 
However, the condition is not sustained in general for cases with interactions 
between spins, and entangled ground states often appear. In the later dis- 
cussion, we are interested in spin-chain models with entangled ground states. 
Let us consider an ordinary case of a correlation function (g\T n O m \g) , not 
decomposed into (g\T n \g)(g\O m \g): 

(g\T n O m \g)^(g\T n \g)(g\O m \g) (10) 

for a certain site n and local operator O m at site m with |n — m| > L + 1. 
Clearly, the state \g) is entangled because all separable ground states satisfy 
(g\T n O m \g) = (g\T n \g) (g\O m \g) . It is proven easily from eq. (|TUJ) that the 
state \g) is not an eigenstate of T n . The reason is as follows. If T n \g) = c n \g) 
is satisfied for a certain real constant c n , the correlation function must be 
written as 

(g\T n O m \g) = c n (g\O m \g) = (g\T n \g)(g\O m \g), 
which contradicts eq. (fTOl) . Therefore we obtain the relation: 

T n \g)^c n \g) (11) 

for arbitrary real constant c n . eq. fill I) gives important information about 
the emergence of a negative energy density as follows. Because the operator 
T n is a Hermitian operator acting on the total Hilbert space of the spin chain, 
T n can be spectrally decomposed into 

T n = ^ e v (n) 1 e v (n) ,k v ,n) {e v (n),k v ,n\, 

where e v {n) are eigenvalues of T n , \e u (n), k Ul n) are corresponding eigenstates 
and the index k u denotes the degeneracy freedom of the eigenvalue e u (n). 
Because {| e u (n), k u , n)} is a complete orthogonal set of basis state vectors of 
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the total Hilbert space, the ground state in the total Hilbert space can be 
uniquely expanded as 

\g) = ^2gu,kA n )\^(n),K,n}. (12) 

u,k v 

Using this expansion, eq. (JTj) gives 

(g\T n \g) = J2^(n)\9,,kM\ 2 = 0. (13) 

v,k v 

Clearly, eq. (|T3|) has no solutions if the lowest eigenvalue e_(n) of T n is 
positive. For the case of e_(n) = 0, eq. (TT5|> has a solution with g , _ ) jj._ (n) 7^ 0. 
The other components g v ^ v ip) vanish. However, this solution implies that 
T n \g) = c n \g) with c n = e_(n) because of eq. ffT2l and apparently contradicts 
eq. ffTTj) . Therefore e_(n) does not take a zero value. Hence, it is concluded 
that e_(n) must be negative: 

e_(n) = - |e_(n)| < 0. 

The average energy density for the corresponding eigenstate |e_(n), k-, n) 
also becomes negative. It is thereby verified that there exist quantum states 
with negative energy density for spin chains satisfying eq. ffTUl) . It should be 
stressed again that even if a state has negative energy density over a certain 
region, there exists positive energy density at other regions and the total 
energy is not negative because of the non-negativity of H. 



3 Local Energy Conservation 



In this section, the local energy conservation of the spin chain system S is 
explained. First, a free system, which is not coupled with external systems, is 
considered. The system evolves by the Hamiltonian in eq. (Jl|). A connected 
site region is introduced, V, given by [n*, nf\ that satisfies nf — rij > 2L — 1. 
The energy of V is defined by 

n } 
n=rn 
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In a free system, the Heisenberg operator of a Schrodinger operator O is 
defined by 0^ H \t) = e ltH Oe~ ltH . The Heisenberg operator of Hy evolves by 
the following equation: 

|<> = i[ir,<>]. (15) 

By substituting a relation given by 

m-l n f +2L 

h= E E + E E T ^ 

n<rii— 2L — 1 n=rii~2L n=rif+l n>rif+2L+l 

into eq. ( fT5l) . it is possible to obtain 

n,-l n/+2L 

(it 



U^ = i { e Tf + e ^.^i- 



n=rn—2L n=ny+l 

Substituting eq. (Tbfj) into the above equation yields the following energy 
conservation relation for an arbitrary initial state 



dt 



WW = J nt -i-J nr 



where the above energy fluxes are given by 

n n+2L 

j n = m\ e T ^ E r i?M ( 16 ) 

m=n— 2L+1 m'=n+l 

If J ni -i = J n , = 0, the energy of V does not change at all. 

Next, a situation is considered in which the spin chain S is locally coupled 
with external systems D and C. S interacts with the external systems only 
at site n Q . The coupling can be switched on and off effectively by the time 
evolution of the switch system C. System D locally controls the spin chain in 
the switch-on interval. The free part of the total Hamiltonian is then written 
as 

H = H + He + Hp, 
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where He is the free Hamiltonian of C, and Hp is of D. In order to capture 
the essence of physics, example cases are considered in which the switch sys- 
tem C is a one-dimensional bosonic Schrodinger field fy(x) with Hamiltonian 



iv f°° 

H c = -- [&(x)d x *(x)-d x &(x)*(xj\dx 

^ J-oo 

1 f°° 

+ 2Mj d ^( x ) d ^( x ) dx 

where v and M are positive real paramters and the operators ty(x) and \&t(x) 
satisfy 

&(x')] = i6(x-x'), 
[V(x), V(x')} = 0, 
[&(x), &(x')] = 0. 

In later discussion, we take the mass parameter M very large. Then v bo- 
comes the propagating velocity of right-moving excitation of the field outside 
the interaction region: 

e itH °^(x)e- itH ° « ^(x - vt). 

The vacuum state |0) of ^ is defined by \l/(x)|0) = 0. By use of the explicit 
example of C, a model of local operation for the spin S at site n Q with energy 
conservation can be made. The total Hamiltonian is given by 



TT TT I Tj( n o) 

n-tot — ti Q + n int , 



Ht ] = gKg> D f &(x)*(x)dx (17) 
J-d 



Here g is a real coupling constant, d fixes the interaction region of \1/ to S 
and D. K^'l is a Hermitian operator acting on the Hilbert space of the 
composite system S and D, and can be decomposed as 

K^; ] D = i no o D + a no -d D , 

— * 

where Op and Op are Hermitian operators acting on the Hilbert space of D. 
In this model, energy conservation is satisfied via 
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d 
dt 



H(t) + H c (t) + H D (t) + H&\t) 



) 







(18) 



where the hat operators denote the Heisenberg operator corresponding to 
Schrodinger operators: 



Switching is realized by scattering a wave packet of \I> via the interaction 
in eq. ( JTTj) . Let us set the initial state of C to a coherent state 



where 4>(x) is amplitude of the coherent state and localized in the space 
region [xo,xi] such that x\ < —d. Let us assume that the total initial state 
is a product state given by |J) = \ip s ) \<j>c), where \ip s ) is the initial 
state of S and \fpjj) is the initial state of D. It is noted that in the initial 
phase, the interaction in eq. (IT7j) is switched off because no wave of ^ exists 
inside the interaction region [—d, d] and H^f* has no contribution. The wave 
packet evolves freely with velocity v until a part of the packet reaches [—d, d]. 
When a part of the wave packet stays in [— d, d], the interaction in eq. ifTTj) 
is effectively switched on and energy can be exchanged among S, C and D. 
In the interval, the wave packet of \& is scattered by the interaction. After 
the scattering process, the scattered wave of \I/ escapes with velocity v from 
the interaction region [—d, d] and H^f 1 again has no contribution. Hence, 
taking account of the disappearance of H^f 1 contribution in the past and 
future of the scattering events, energy conservation in eq. (lisp clearly yields 
the following relation: 



where U is the initial time, and tf is the final time of the scattering process. 
It should be stressed that the scattering interval tf — U can go to zero as v 
and g approach infinity. It should be recalled that the nonrelativistic limit is 
assumed in this paper, and the speed of light is assumed to be infinity. Thus, 
v can also approach infinity. At this limit, the time evolution of the spins of 
S at all sites except site n Q can be neglected during the short-time scattering 



0(t) = e itHtot Oe 




(I\ (H(U) + H C {U) + Hd{U)) |/> = </| (#(*/) + H c (tf) + H D {t f )) \I) 



(19) 
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process. Then, energy exchange occurs not globally but locally among spins 
around site n Q , C and D. Therefore, eq. fll9p can be rewritten as 

n +L \ 

f n (U)+H c (U) + H D (U)\ \I) 

=n —L / 
n +L \ 

(J|l T n (tf)+Hc(tf)+Hn(tf)\\I). (20) 

\n=n -L / 

If (I\ (H(tf) — H(ti) J \I) is positive, the spin chain gets energy from C and 



D in the short-time scattering interval. Inversely, if (I\ (jk(tj) — H(ti)J \I) 

is negative, the spin chain gives energy to the composite system of C and D. 
This fact becomes significant in the QET protocol proposed later. 

4 Local Measurement of the Ground State 
and Entanglement Breaking 

For the QET protocol, it is important to consider the local measurement of 
an entangled ground state \g) of a spin chain. Let us consider a spin at site 
riA in the spin chain and a Hermitian unitary local operator o a = "a • &n A - 
Here ua is a three-dimensional real unit vector and a nA is the Pauli spin 
matrix vector at site Ua- The eigenvalues of a a are (— 1) M with fi — 0,1. The 
spectral expansion of a a is given by 

a A — UA-a A — Yl Pa (^) ' ( 21 ) 

m=o,i 

where Pa(h) is a projective operator onto the eigensubspace with \i. Consider 
that Alice at site ua performs a projective measurement of a a for the ground 
state \g) and assume that she obtains the measurement result \i. The post- 
measurement state with \i is given by 

Pa(m)\9), (22) 



where Pa{^) = (pI-Pa^) Therefore, the average post-measurement state 
is a mixed quantum state p' given by 
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p' = pa(h)—f=^ Pa ^) \g){g\PA(n) —r=r^ 
M t^i vp aw vp aw 

= Y. p A^)\g){g\P A {v)- 

M=0,1 

It should be noted that p' is a quantum state which satisfies 

Tr[p'} = Tr[\g)(g\], (23) 

n A n A 

where Tr nA means a partial trace in terms of the measured spin at site ha- 
Therefore, the quantum fluctuation of p' is the same as that of the ground 
state except at site tia- It should also be noted that the above projective 
measurement needs energy. Because the ground state has zero energy, the 
input Ea is calculated as follows: 

E A = Tr [p'H] - (g\H\g) = £ (g\P A (ji) HP A (/i) \g). (24) 

M=0,1 

Because of the non-negativity of H, Ea is non-negative. In general, Pa (p) \g) 
is not proportional to \g) due to entanglement and we obtain a positive value 
of Ea- Because the measurement is performed locally, the excitation energy 
of the spin chain is localized around site ha soon after the measurement. For 
convenience, we introduce a localized energy operator around site n by 

n+L 
m=n—L 

H n includes all T m 's for which expectation values change if quantum oper- 
ations act on site n. From eq. (125|) . it can be shown that 

Tr [p'H nA ] = Ea- 

Also, it is proven that 

Tr [p'T n } = (26) 

where |n — ua\ > L. eq. f[2"6"|) implies localization of the energy distribution 
around site tia- 
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Next, consider that Alice at site ua attempts to completely withdraw the 
input energy E\ by local operations soon after the measurement. Her avail- 
able processes are expressed by arbitrary trace-preserving completely-positive 
(TPCP) maps Y a of the spin state at site ua- In spite of any use of local 
TPCP maps, her attempt cannot be achieved because the local measurement 
breaks the entanglement between the spin at site ua and different-site spins 
in the spin chain. If Alice wants to recover the original ground state, she 
must recreate the entanglement broken by the first measurement. However, 
entanglement generation needs nonlocal operations in general. Therefore, 
she cannot recover the ground state perfectly by her local operations alone. 
This observation implies that there is nonvanishing residual energy E res in 
Alice's local cooling process to extract the input energy: 

E res = mm(Tr[r A [p']H nA })>0. 

Therefore, interestingly, Alice is not able to use the residual energy even 
though it is right in front of her. In the next section, we argue that a part 
of the residual energy can be extracted by Bob, who is situated apart from 
Alice and performs local operations dependent on the measurement result \i 
of a a to a spin on the chain in front of him. This protocol is QET. 



5 Quantum Energy Teleportation 



In this section, we discuss QET for spin chain systems. As mentioned in 
the introduction, we concentrate on short time scales, in which dynamical 
evolution induced by the Hamiltonian of the spin chain is negligible. On the 
other hand, we also assume the nonrelativistic limit that LOCC for the spins 
can be repeated many times even in a short time interval. Initially, Alice and 
Bob share many near-critical spin chains in the ground state \g), which is 
entangled and has a large correlation length /. Alice is situated at site ua and 
Bob at site ns- Alice is a good distance from Bob: \n A — n B \ ~ 0(1) 3> 1. 
Consider a protocol consisting of the following steps. (I) Alice performs a 
measurement of a local observable of the spin at site n A and obtains the 
measurement result \i. In the measurement process, the spin chain is locally 
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excited with some energy input. (II) Alice announces to Bob the result /x by 
a classical channel at the speed of light or near the light velocity. Because the 
time interval of the annoucement is very short, the system does not evolve 
and no energy flow appears between Alice and Bob. (Ill) Bob performs 
a local operation dependent on /x to the spin at site %. QET is defined 
by this protocol if the expectation value of H nB is negative after the local 
operation of Bob. Because the expectation value of H ng is zero just before 
Bob's operation and local energy conservation holds around Bob's region as 
mentioned in section 3, the negative value after Bob's operation implies that 
energy release occurs around site Ub to external systems during the operation 
interval. (For the case in section 3, the energy moves from the spin chain S 
to the external systems C and D.) If Alice inputs no energy into the spin 
chain or Bob does not use the information about /x, Bob is not able to extract 
energy from the spin chain at all. Therefore, it can be said that a part of 
the energy input by Alice's measurement at site ua is effectively transported 
to site «b by LOCC. If we select as the initial state of the spin chain, not 
the ground state but an excited state which has nonzero energy distribution 
around Bob, it is not unusual that Bob's local operation can extract energy 
from the spin chain. However, in QET, the ground state is actually selected 
as the initial state and the local energy extraction is nontrivial. 

Next we propose an explicit protocol of QET. The protocol with L = 1 
is illustrated in Figure 1. 

(I) Alice performs a projective measurement of the observable a a in eq. 
(121]) for the ground state \g) and obtains the measurement result /x = or 
1. Alice inputs energy Ea in eq. (I24p to the spin chain system in order to 
achieve the local measurement, as seen in section 4. 

(II) Alice announces to Bob the result /x by a classical channel. In the 
announcement process, both time evolution of the system and emergence of 
energy flux do not happen because the time interval is assumed very short. 

(III) To a spin at site tib, Bob performs a local unitary operation Vb (/x) 
depending on the value of /x, which is defined by 

V B (fi) = 1 cos 9 + i(-lfa B sin9. (27) 
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Here <Jb = ub ■ o nB , ub is a three-dimensional real unit vector and a ng is 

the Pauli spin matrix vector at site n#. The real parameter 9 is defined by 

cos (29) = 7 =JL= , (28) 

sin(2£) = ~^==. (29) 

where the real parameters £ and 77 are given by 

(=(g\<j B H<TB\g), (30) 

v = (gWA^B\g) (3i) 

with 

cr B = z [#„ B , a B ] ■ (32) 
For the case with 77 7^ 0, Bob obtains positive energy: 

' ' (33) 



on average from the spin chain in the process of the local operation. 

The value of £ is non-negative due to the non-negativity of H. Because 
of the hermicity and commutativity of a a and &b, the reality of i] is easily 
proven as follows: 

v* = (g\vA<7B\g)* = {g\(? B a A \g) 
= (g\(r A & B \g} = 77. 

Using [T n , ctb] = for |n — ub\ > L, eq. fl32l) can be rewritten as 

b B = i [H, a B ] , (34) 

and thus a b can be interpreted as a time- derivative operator of the Heisen- 
berg operator exp [itH] Ob exp[— itH] at t = 0. After Bob's operation, the 
post-measurement state in eq. (1221) is transformed into 

-^=VB^)PA^)\g). 
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Here we have neglected time evolution of the spin chain between Alice's 
measurement and Bob's operation assuming the relation of eq. (JTJ) or eq. 
(T5]). The average state is given by 



]T p A (p)—L=V B (jjl) P A (p) \g) (g\P A (jjl) V B (p) - 7 =L= 
£ V B (p)P A (p) \g)(g\P A (»)V B (p). 



^=0,1 



It is straightforward to calculate the average localized energy Tr [pH nB ] after 
Bob's operation, as follows. Firstly, by the commutativity of P A (p) and 
V B (p) H nB V B (p), we obtain 

Tr [pH nB ] = Y,^\PA^){v B (p)H nB V B (p))P A (p)\g) 

M=0,1 

= J2<9\ p A(^(vU^H nB V B (fM))\g), 

where we have used P A (p) 2 = P A (p). Substituting eq. ( l27l) into the above 
relation yields the following expression of Tr [pH nB ]: 

Tr [ P H nB ] 

= £ (9\Pa fa) (I cosO -i (-If a B sin 9) H ng (I cos9 + i (-If a B sin 9) \g) 

m=o,i 

= cos 2 %| ( J2 P A(»))Hn B \g)+s™ 2 9(9\ [J2 PA^))cj B H nB a B \g) 
V=o,i / V=o,i / 

+zcos0sin%| ( (-^Pa^)) [Hn B , cr B ]\g). 

\M=0,1 / 

By taking account of the completeness relation of P A (p) and the spectral 
decomposition of a A in eq. ( I2TT) . we get the following relation: 

Tr[pH nB ] = cos 2 9 (g\H nB \g) + sin 2 9 (g\a B H nB cr B \g) 
+i cos 9 sin 9 (g\a A [H nB , a B )\g). 
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Because (g\H nB \g) = is satisfied, the form given by 

Ti [pH nB ] = {g\(y B H nB a B \g) sin 2 9 + r] cos 9 sin 9 (35) 
is obtained. Here, the relation: 

(g\a B H ng a B \g) = (g\a B I H - T n ,\a B \g) 

\ ri£[n B -L, n B +L] J 

= (g\vBHa B \g) 

is satisfied due to the following relation: 

(g\(T B T n ,(T B \g) = (g\cr B (x B T nl \g} = (g\T n >\g) = 

for n' ^ [n B — L, n B + L}. Thus, eq. (!35j) is rewritten as 

Tr [pH nB \ = £ sin 2 6* + 77 cos 6* sin 9 

= |(l-cos(20)) + |sin(20). (36) 

The parameter 9 is now fixed to make Tr [pH nB ] as negative as possible. 
This can be achieved by taking 9 as in eqs. ( 1281) and ( [291) . Tr [p-ffnj is then 
evaluated as 



Tr [pH nB ] = \ 



(37) 



by substituting eqs. (1281) and (|29|) into the relation derived from eq. (!36|) . If 
?7 7^ 0, it is clear that Tr [p-ff nB ] m eq. (1371) is negative: 

Tr [p# nB ] < 0. (38) 

eq. (13"51) is a significant result. Before step (III), the energy around Bob is 
zero: 

Tr [p'H nB \ = 0. 

After Bob's local operations in step (III), the localized energy around site n B 
becomes negative. Respecting local energy conservation explained in section 
3, this means that positive energy E B must be emitted from the spin chain 
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to Bob (more precisely, to Bob's devices to perform Vg (//)). More explicitly, 
local energy conservation around site ub yields the relation: 



E B + Tr [pH nB ] = Tr [p'H nB ] = 0. 



From this relation, we finally obtain eq. (1331) as follows: 



E B 



Tr [p'H nB ] - Tr [pH nB ] 



As seen in eq. ( IHTj) . 77 is given by a two-point correlation function of (semi- 
) local operators. If the ground state is separable, it turns out from eq. flMl) 
and H\g) = that 77 vanishes as follows: 



Therefore, it can easily be checked from eq. that Bob gains no energy, 
as should be the case. However, the correlation function 77 does not van- 
ish in general because entanglement yields correlations among the quantum 
fluctuations at each point. For example, the critical Ising spin chains with 
transversal magnetic field have nonvanishing 77 [6J. In general, any spin chain 
models without fine-tuning may have nonvanishing 77, just like the Ising spin 
chains, and the protocol in this paper becomes effective barring fine-tuned ex- 
ceptions. Schematic plots of the expectation values of energy density for the 
steps of the QET are given in Figs. 2-4. From these figures, it can be easily 
understood that there exists no energy flow between Alice and Bob. Energy 
is locally transported from the spin chain to Bob's device for the operation 
by simultaneously generating negative energy — E B in the spin chain and 
positive energy +Eb in the device of Bob. Hence, local energy conservation 
is exactly maintained in the protocol. 

It is worth commenting here that the ground state \g) has a typical cor- 
relation site length I, over which the correlation between two sites and the 
value of 77 decay rapidly. Hence, the protocol is more effective in teleporting 
energy between Alice and Bob with \ua — Ub\ ~ I, because the energy gain of 
Bob increases when 77 is large. To achieve long-range QET, it is preferable to 
choose near-critical spin chain systems which have entangled ground states 
with long-range correlation, that is, I 3> 1. For example, the near-critical 



n 



(g\<7A0B\g) = (gWA\g)(g\vB\g) 

i(g\a A \g)(g\(Ha B - <r B H)\g) = 0. 
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Ising spin chains in the presence of transversal magnetic field with b ~ h > 
in eq. (jSD may be good candidates. 

In the above analysis, we have shown that Bob obtains energy from the 
spin chain. However, even after the last step of the protocol, there exists 
energy Ea that Alice had to first deposit to the spin chain by her measure- 
ment. Actually, using H nj y B (p) = V B (p) H nA , we are able to check easily 
that Tr [pH nA ] = Tr [p'H nA \. Let us imagine that Alice attempts to com- 
pletely withdraw Ea by local operations after step (III). If this was possible, 
the energy gain Eb of Bob might have no cost. However, if so, the total en- 
ergy of the spin chain becomes equal to — Eb and negative. Meanwhile, we 
know that the total energy of the spin chains must be nonnegative. Hence, 
Alice cannot withdraw energy larger than Ea — E B by local operations at 
site ha- The reason of existence of the residual energy has been already ex- 
plained in section 4 from the viewpoint of entanglement breaking. It may be 
instructive to stress the simple fact that if the deposited energy Ea vanishes, 
Eb also vanishes. This implies that Bob releases a part of energy poured by 
Alice's measurement. 

6 Conclusion 

In this paper, a protocol of QET is proposed in spin chain systems that 
respects all physical laws including causality and local energy conservation. 
Energy is effectively transported from Alice to Bob by only LOCC, using 
local excitations with negative energy and ground-state entanglement. The 
energy input by Alice is given by eq. (I24I) and the energy gain by Bob is 
given by eq. (|33|) . 

If we do not consider the ground state, but rather an excited state of 
the spin chain which has nonzero energy distribution around Bob, it is not 
anomalous that Bob can extract energy from the spin chain. In the QET 
protocol, the initial state is set to be the ground state and energy extraction 
by Bob is nontrivial. At first glance, energy extraction from the ground state 
by Bob appears impossible due to local energy conservation because Alice 
only excites locally the spin chain by her measurements. Based on this, we 
may wonder why energy is teleported in QET without any physical agent in 
the protocol. The answer may be briefly summarized as follows. As stressed 
in section 2, quantum mechanics allows negative-energy- density regions by 
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controlling quantum fluctuation. Hence, even if Bob has no energy on aver- 
age around him, the value of the energy density around him can decrease. 
Therefore, the ground state can be locally regarded as "an excited state", 
compared with the state with negative energy around Bob. In the QET pro- 
tocol, Bob extracts this local excitation energy by using Alice's measurement 
result. In this way, it can be said that the energy, which Bob will obtain, 
existed around Bob before the start of the QET protocol. Therefore, we do 
not need any transfer of energy from Alice to Bob. Of course, the energy 
hidden behind Bob's region is not always available. In the protocol, classical 
information about the measurement result of Alice becomes a key to extract 
energy by Bob. Without knowledge about the measurement result, Bob can- 
not get any energy from the local ground state. The essential reason why the 
measurement result allows Bob to get energy is as follows. Because the spin 
of Alice is entangled with spins around Bob, the measurement result includes 
information about the quantum fluctuations around Bob. Bob infers from 
the measurement information how the quantum fluctuations behave around 
him and can choose a unitary operation to extract but not give energy from 
the spin chain. Consequently, the protocol is able to transport energy by 
LOCC without breaking local energy conservation and without the existence 
of any excited physical entity. 
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Figure Caption 



Figure 1: QET protocol 

(1) The eigenvalues of ua • &a are (— 1) M with fj, — 0,1. Alice performs 
a local projective measurement of the observable u A ■ a A to her spin in the 
ground state \g) and obtains the measurement result /i. Alice must input 
energy Ea to the spin chain system in order to achieve the local measurement. 
(2) Alice announces to Bob the result \x by a classical channel. (3) Bob 
performs a local unitary operation Vb (fj) to his spin, depending on the value 
of /i. Bob obtains energy output E B on average from the spin chain system 
in the process of the local operation. 

Figure 2: Schematic plot of the expectation value of energy density in 
Step (I). Alice measures ua • <$a and obtains the result \x. The energy Ea is 
poured into the spin chain by the local measurement. 

Figure 3: Schematic plot of the expectation value of energy density in 
Step (II). Alice announces fj, to Bob. Bob perfomes result-dependent unitary 
operation to his spin. 

Figure 4: Schematic plot of the expectation value of energy density in Step 
(III). Negative energy — Eb appears in the spin chain and positive energy 
+Eb is released outside. Bob is able to use the energy. 
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